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Motivated by the heavy ion collision experiments there is much activity in studying the hydrodynamical properties of 
non-abelian (quark-gluon) plasma's. A major question is how to deal with color currents. Although not widely appreciated, 
quite similar issues arise in condensed matter physics in the context of the transport of spins in the presence of spin-orbit 
coupling. The key insight is that the Pauli Hamiltonian governing the leading relativistic corrections in condensed matter 
systems can be rewritten in a language of SU (2) covariant derivatives where the role of the non-abelian gauge fields is taken 
by the physical electromagnetic fields: the Pauli system can be viewed as Yang-Mills quantum mechanics in a 'fixed frame', 
and it can be viewed as an 'analogous system' for non abelian transport in the same spirit as Volovik's identification of the 
He superfiuids as analogies for quantum fields in curved space time. We take a similar perspective as Jackiw and coworkers in 
their recent study of non-abelian hydrodynamics, twisting the interpretation into the 'fixed frame' context, to find out what 
this means for spin transport in condensed matter systems. We present an extension of Jackiw's scheme: non-abelian 
hydrodynamical currents can be factored in a 'non-coherent' classical part, and a coherent part requiring macroscopic 
non-abelian quantum entanglement. Hereby it becomes particularly manifest that non-abelian fiuid fiow is a much richer 
afi'air than familiar hydrodynamics, and this permits us to classify the various spin transport phenomena in condensed matter 
physics in an unifying framework. The "particle based hydrodynamics" of Jackiw et al. is recognized as the high temperature 
spin transport associated with semiconductor spintronics. In this context the absence of faithful hydrodynamics is well 
known, but in our formulation it is directly associated with the fact that the covariant conservation of non-abelian currents 
turns into a disastrous non-conservation of the incoherent spin currents of the high temperature limit. We analyze the 
quantum-mechanical single particle currents of relevance to mesoscopic transport with as highlight the Ahronov-Casher effect, 

where we demonstrate that the intricacies of the non-abelian transport render this effect to be much more fragile than its 
abelian analogue, the Ahronov-Bohm effect. We subsequently focus on spin flows protected by order parameters. At present 
there is much interest in multiferroics where non-coUinear magnetic order triggers macroscopic electric polarization via the 
spin-orbit coupling. We identify this to be a peculiarity of coherent non-abelian hydrodynamics: although there is no net 
particle transport, the spin entanglement is transported in these magnets and the coherent spin 'super' current in turn 
translates into electric flelds with the bonus that due to the requirement of single valuedness of the magnetic order parameter 
a true hydrodynamics is restored. Finally, 'flxed-frame' coherent non-abelian transport comes to its full glory in spin-orbit 
coupled 'spin superffuids', and we demonstrate a new effect: the trapping of electrical line charge being a fixed frame, 
non-abelian analogue of the familiar magnetic fiux trapping by normal superconductors. The only known physical examples of 
such spin superfiuids are the ''He A- and B phase where unfortunately the spin-orbit coupling is so weak that it appears 

impossible to observe these effects. 



PACS numbers: 73.43.-f,72.25.Dc,72.25.Hg 

I. INTRODUCTION 

It is a remarkable development that in various branches 
of physics there is a revival going on of the long standing 
problem of how non- Abelian entities are transported over 
macroscopic distances. An important stage is condensed 
matter physics. A first major development is spintron- 
ics, the pursuit to use the electron spin instead of its 
charge for switching purposesii^ii^ii^i^, with a main focus 
on transport in conventional semiconductors. Spin-orbit 
coupling is needed to create and manipulate these spin 
currents , and it has become increasingly clear that trans- 
port phenomena are possible that are quite different from 
straightforward electrical transport. A typical example is 
the spin-Hall effec1jii2,i^, defined through the macroscopic 
transport equation, 

JT = OsH^ialEl (1) 

where eiai is the 3-dimensional Levi-Civita tensor and Ei 
is the electrical field. The specialty is that since both 
j'f and El are even under time reversal, the transport 
coefficient asH is also even under time reversal, indicat- 
ing that this corresponds with a dissipationless transport 



phenomenon. An older development is the mesoscopic 
spin-transport analogue of the Aharonov-Bohm effect, 
called the Aharonov-Casher effect^: upon transversing 
a loop containing an electrically charged wire the spin 
conductance will show oscillations with a period set by 
the strength of the spin-orbit coupling and the enclosed 
electrical line-charge. 

A rather independent development in condensed mat- 
ter physics is the recent focus on the multiferroics. This 
refers to substances that show simultaneous ferroelectric- 
and ferromagnetic order at low temperatures, and these 
two different types of order do rather strongly depend 
on each other. It became clear recently that at least in 
an important subclass of these systems one can explain 
the phenomenon in a language invoking dissipationless 
spin transport^i^: one needs a magnetic order charac- 
terized by spirals such that 'automatically' spin currents 
are flowing, that in turn via spin-orbit coupling induce 
electrical fields responsible for the ferroelectricity. 

The final condensed matter example is one that was 
lying dormant over the last years: the superfiuids real- 
ized in ^He. A way to conceptualize the intricate order 
parameters of the A- and B-phase^°'^^ is to view these as 
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non-Abelian ('spin-like') superfluids. The intricacies of 
the topological defects in these phases is of course very 
well known, but matters get even more interesting when 
considering the effects on the superflow of macroscopic 
electrical fields, mediated by the very small but finite 
spin-orbit coupling. This subject has been barely stud- 
ied: there is just one paper by Mineev and Volovik^^ 
addressing these matters systematically. 

A very different pursuit is the investigation of the 
quark-gluon plasma's presumably generated at the 
Brookhaven heavy-ion collider. This might surprise the 
reader: what is the relationship between the flow of spin 
in the presence of spin-orbit coupling in the cold con- 
densed matter systems and this high temperature QCD 
affair? There is actually a very deep connection that 
was already realized quite some time ago. Goldhaber— 
and later Froehlich et a?"*, Balatskii and Altshuler— and 
others realized that in the presence of spin-orbit cou- 
pling spin is subjected to a parallel transport principle 
that is quite similar to the parallel transport of matter 
fields in Yang-Mills non-Abelian gauge theory, underly- 
ing for instance QCD. This follows from a simple rewrit- 
ing of the Pauli-equation, the Schroedinger equation tak- 
ing into account the leading relativistic corrections: the 
spin-fields are just subjected to covariant derivatives of 
the Yang-Mills kind, see Eq.'s ([S]),®. However, the dif- 
ference is that the 'gauge' fields appearing in these co- 
variant derivatives are actually physical fields. These are 
just proportional to the electrical- and magnetic fields. 
Surely, this renders the problem of spin transport in con- 
densed matter systems to be dynamically very different 
from the fundamental Yang-Mills theory of the standard 
model. However, the parallel transport structure has a 
'life of its own': it implies certain generalities that are 
even independent of the 'gauge' field being real gauge or 
physical. 

For all the examples we alluded to in the above, 
one is dealing with macroscopic numbers of particles 
that are collectively transporting non-Abelian quantum 
numbers over macroscopic distances and times. In the 
Abelian realms of electrical charge or mass a univer- 
sal description of this transport is available in the form 
of hydrodynamics, be it the hydrodynamics of water, 
the magneto-hydrodynamics of charged plasma's, or the 
quantum-hydrodynamics of superfiuids and supercon- 
ductors. Henceforth, to get anywhere in terms of a sys- 
tematic description one would like to know how to think 
in a hydrodynamical fashion about the macroscopic fiow 
of non-Abelian entities, including spin. 

In the condensed matter context one finds pragmatic, 
case to case approaches that are not necessarily wrong, 
but are less revealing regarding the underlying 'universal' 
structure: in spintronics one solves Boltzmann transport 
equations, limited to dilute and weakly interacting sys- 
tems. In the quark-gluon plasma's one find a similar 
attitude, augmented by RPA-type considerations to deal 
with the dynamics of the gauge fields. In the multifer- 
roics one rests on a rather complete understanding of the 



order parameter structure. 

The question remains: what is non-Abelian hydro- 
dynamics? To the best of our knowledge this issue is 
only addressed on the fundamental level by Jackiw and 
coworkersi^iii and their work forms a main inspiration 
for this review. The unsettling answer seems to be: non- 
Abelian hydrodynamics in the conventional sense of de- 
scribing the collective flow of quantum numbers in the 
classical liquid does not even exist! The impossibility to 
define 'soft' hydrodynamical degrees of freedom is rooted 
in the non-Abelian parallel transport structure per se and 
is therefore shared by high temperature QCD and spin- 
tronics. 

The root of the trouble is that non-Abelian currents 
do not obey a continuity equation but are instead only 
covariantly conserved, as we will explain in detail in sec- 
tion |Vl It is well known that covariant conservation laws 
do not lead to global conservation laws, and the lack of 
globally conserved quantities makes it impossible to deal 
with matters in terms of a universal hydrodynamical de- 
scription. This appears to be a most serious problem for 
the description of the 'non-Abelian fire balls' created in 
Brookhaven. In the spintronics context it is well known 
under the denominator of 'spin relaxation': when a spin 
current is created, it will plainly disappear after some 
characteristic spin relaxation determined mostly by the 
characteristic spin-orbit coupling strength of the mate- 
rial. 

In this review we will approach the subject of spin 
transport in the presence of spin-orbit coupling from the 
perspective of the non-Abelian parallel transport princi- 
ple. At least to our perception, this makes it possible to 
address matters in a rather unifying, systematical way. 
It is not a-priori clear how the various spin transport 
phenomena identified in condensed matter relate to each 
other and we hope to convince the reader that they are 
different sides of the same non-Abelian hydrodynamical 
coin. Except for the inspiration we have found in the pa- 
pers by Jackiw and coworker o^^'^^ we will largely ignore 
the subject of the fundamental non-Abelian plasma, al- 
though we do hope that the 'analogous systems' we iden- 
tify in the condensed matter system might form a source 
of inspiration for those working on the fundamental side. 

Besides bringing some order to the subject, in the 
course of the development we found quite a number of 
new and original results that are consequential for the 
general, unified understanding. We will start out on the 
pedestrian level of quantum- mechanics (section III) , dis- 
cussing in detail how the probability densities of non- 
Abelian quantum numbers are transported by isolated 
quantum particles and how this relates to spin-orbit cou- 
pling (Section IV). We will derive here equations that 
are governing the mesoscopics, like the Aharonov-Casher 
(AC) effect, in a completely general form. A main con- 
clusion will be that already on this level the troubles with 
the macroscopic hydrodynamics are shimmering through: 
the AC effect is more fragile than the Abelian Aharonov- 
Bohm effect, in the sense that the experimentalists have 
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to be much more careful in designing their machines in 
order to find the AC signal. 

In the short section V we revisit the non-Abelian co- 
variant conservation laws, introducing a parametriza- 
tion that we perceive as very useful: different from the 
Abelian case, non-Abelian currents can be viewed as be- 
ing composed of both a coherent, 'spin' entangled part 
and a factorisable incoherent part. This difference is at 
the core of our classification of non-Abelian fluids. The 
non-coherent current is responsible for the transport in 
the high temperature liquid. The coherent current is re- 
sponsible for the multiferroic effects, the Meissner 'dia- 
magnetic' screening currents in the fundamental non- 
Abelian Higgs phase, but also for the non-Abelian su- 
percurrents in true spin-supcrfluids like the '^He A- and 
B phase. 

The next step is to deduce the macroscopic hydrody- 
namics from the microscopic constituent equations and 
here we follow Jackiw et. ali^'^ closely. Their 'parti- 
cle based' non-Abelian hydrodynamics is just associated 
with the classical hydrodynamics of the high temperature 
spin-fluid and here the lack of hydrodynamical descrip- 
tion hits full force: we hope that the high energy physi- 
cists flnd our simple 'spintronics' examples illuminating 
(Section VI). 

After a short technical section devoted to the workings 
of electrodynamics in the SO problem (section VII), we 
turn to the 'super' spin currents of the multiferroics (Sec- 
tion VIII) . As we will show, these are rooted in the coher- 
ent non-Abelian currents and this renders it to be quite 
similar but subtly different from the 'true' supercurrents 
of the spin superfluid: it turns out that in contrast to 
the latter they can create electrical charge! This is also 
a most elementary context to introduce a notion that we 
perceive as the most important feature of non-Abelian 
fluid theory. In Abelian hydrodynamics it is well under- 
stood when the superfluid order sets in, its rigidity does 
change the hydrodynamics: it renders the hydrodynam- 
ics of the superfluid to be irrotational having the twofold 
effect that the circulation in the superfluid can only occur 
in the form of massive, quantized vorticity while at low 
energy the superfluid is irrotational so that it behaves like 
a dissipationless ideal Euler liquid. In the non-Abelian 
fluid the impact of the order parameter is more dramatic: 
its rigidity removes the multivaluedness associated with 
the covariant derivatives and hydrodynamics is restored! 

This bring us to our last subject where we have most 
original results to offer: the hydrodynamics of spin- 
orbit coupled spin-superfluids (Section IX) . These are the 
'flxed frame' analogs of the non-Abelian Higgs phase and 
we perceive them as the most beautiful physical species 
one encounters in the non-Abelian fluid context. Un- 
fortunately, they do not seem to be prolific in nature. 
The "^He superfluids belong to this category but it is 
an unfortunate circumstance that the spin-orbit coupling 
is so weak that one encounters insurmountable difficul- 
ties in the experimental study of its effects. Still we will 
use them as an exercise ground to demonstrate how one 



should deal with more complicated non-Abelian struc- 
tures (Section XI), and we will also address the issue of 
where to look for other spin-superfluids in the concluding 
section (Section XII). 

To raise the appetite of the reader let us start out 
presenting some wizardry that should be possible to re- 
alize in a laboratory when a spin-superfiuid would be dis- 
covered with a sizable spin-orbit coupling: how the elu- 
sive spin-superfluid manages to trap electrical line charge 
(section II), to be explained in detail in Section X. 



II. THE APPETIZER: TRAPPING QUANTIZED 
ELECTRICITY. 

Imagine a cylindrical vessel, made out of plastic while 
its walls are coated with a thin layer of gold. Through 
the center this vessel a gold wire is threaded and care is 
taken that it is not in contact with the gold on the walls. 
Fill this container to the brim with a putative liquid that 
can become a spin superfluid (liquid ^He would work if 
it did not contain a dipolar interaction that voids the 
physics ) in its normal state and apply now a large bias 
to the wire keeping the walls grounded, see Fig. [TJ Since 
it is a capacitor, the wire will charge up relative to the 
walls. Take care that the line charge density on the wire 
is pretty close to a formidable 2.6 x 10~^ Coulomb per 
meter in the case that this fluid would be like ^He . 
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FIG. 1: A superfluid ^He container acts as a capacitor capa- 
ble of trapping a quantized electrical line charge density via 
the electric field generated by persistent spin Hall currents. 
This is te analog of magnetic flux trapping in superconductors 
by persistent charge supercurrents. 

Having this accomplished, cool the liquid through its 
spin superfluid phase transition temperature Tc- Remove 
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now the voltage and hold the end of the wire close to 
the vessel's wall. Given that the charge on the wire is 
huge, one anticipates a disastrous decharging spark but 
.... nothing happens! 

It is now time to switch off the dilution fridge. Upon 
monitoring the rising temperature, right at Tc where the 
spin superfluid turns normal a spark jumps from the wire 
to the vessel, grilling the machinery into a pile of black 
rubble. 

This is actually a joke. In Section X we will present 
the theoretical proof that this experiment can actually 
be done. There is a caveat, however. The only substance 
that has been identified, capable of doing this trick is he- 
lium III were it not for the dipolar interaction preventing 
it being the desired spin superfluid. But even if we were 
God and we could turn the dipolar locking to zero mak- 
ing Helium III into the right spin superfluid, there would 
still be trouble. In order to prevent bad things to hap- 
pen one needs a vessel with a cross sectional area that is 
roughly equal to the area of Alaska. Given that there is 
only some 170 kg of helium on our planet, it occurs that 
this experiment cannot be practically accomplished. 

What is going on here? This effect is analogous to 
magnetic flux trapping by superconducting rings. One 
starts out there with the ring in the normal state, in the 
presence of an external magnetic field. One cycles the 
ring below the transition temperature, and after switch- 
ing off the external magnetic field a quantized magnetic 
flux is trapped by the ring. Upon cycling back to the 
normal state this flux is expelled. Read for the magnetic 
flux the electrical line charge, and for the electrical super- 
conductor the spin-supcrfluid and the analogy is clear. 

This reveals that in both cases a similar parallel trans- 
port principle is at work. It is surely not so that this can 
be understood by simple electro-magnetic duality: the 
analogy is imprecise because of the fact that the physical 
field enters in the spin-superfluid problem via the spin- 
orbit coupling in the same way the vector potential enters 
in superconductivity. This has the ramification that the 
electrical monopole density takes the role of the magnetic 
flux, where the former takes the role of physical incarna- 
tion of the pure gauge Dirac string associated with the 
latter. 

The readers familiar with the Aharonov-Casher effect 
should hear a bell ringing^^. This can indeed be con- 
sidered as just the 'rigid' version of the AC effect, in 
the same way that flux trapping is the rigid counterpart 
of the mesoscopic Aharonov-Bohm effect. On the sin- 
gle particle level, the external electromagnetic fields pre- 
scribe the behavior of the particles, while in the ordered 
state the order parameter has the power to impose its 
will on the electromagnetic fields. 

This electrical line-charge trapping effect summarizes 
neatly the deep but incomplete relations between real 
gauge theory and the working of spin-orbit coupling. It 
will be explained in great detail in sections IX and X, but 
before we get there we first have to cross some terrain. 



III. QUANTUM MECHANICS OF SPIN-ORBIT 
COUPLED SYSTEMS 



To address the transport of spin in the presence of 
spin-orbit (SO) coupling we will follow a strategy well 
known from conventional quantum mechanical transport 
theory. We will first analyze the single particle quantum- 
mechanical probability currents and densities. The start- 
ing point is the Pauli equation, the generalization of the 
Schrodinger equation containing the leading relativistic 
corrections as derived by expanding the Dirac equation 
using the inverse electron rest mass as expansion param- 
eter. We will first review the discovery by Volovik and 
Minee^ii^, Balatskii and Altshuler— and Froehlich and 
othersi^ of the non-Abelian parallel transport structure 
hidden in this equation, to subsequently analyze in some 
detail the equations governing the spin-probability cur- 
rents. In fact, this is closely related to the transport of 
color currents in real Yang-Mills theory: the fact that in 
the SO problem the 'gauge fields' are physical fields is 
of secondary importance since the most pressing issues 
regarding non-Abelian transport theory hang together 
with parallel transport. For these purposes, the spin- 
orbit 'fixed-frame' incarnation has roughly the status as 
a representative gauge fix. In fact, the development in 
this section has a substantial overlap with the work of 
Jackiw and co-workers dedicated to the development of 
a description of non-Abelian fluid dynamic o^^i^^ . We per- 
ceive the application to the specific context of SO coupled 
spin fluid dynamics as clarifying and demystifying in sev- 
eral regards. We will identify their 'particle based' fluid 
dynamics with the high temperature, classical spin fluid 
where the lack of true hydrodynamics is well established, 
also experimentally. Their 'field based' hydrodynamics 
can be directly associated with the coherent superflows 
associated with the SO coupled spin superfluids where at 
least in equilibrium a sense of a protected hydrodynami- 
cal sector is restored. 

The development in this section have a direct relevance 
to mesoscopic transport phenomena (like the Aharonov- 
Casher effects^' but here our primairy aim is to set 
up the system of microscopic, constituent equations to 
be used in the subsequent sections to derive the various 
macroscopic fluid theories. The starting point is the well 
known Pauli-equation describing mildly relativistic par- 
ticles. This can be written in the form of a Lagrangian 
density in terms of spinors ip, 
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where 



E = -VAo - doA , 



S = V X A 
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are the usual U{1) gauge fields associated with the 
electromagnetic fields E and B. The relativistic correc- 
tions are present in the terms containing the quantity 
q, proportional to the Bohr magneton, and the time-like 
first term cx i? is the usual Zeeman term while the space- 
like terms oc E corresponds with spin-orbital coupling. 

The recognition that this has much to do with a non- 
Abelian parallel transport structure, due to Mineev and 
Volovik^^, Goldhaber— and Froehlich et al.^^ is in fact 
very simple. Just redefine the magnetic- and electric field 
strengths as follows. 
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Define covariant derivatives as usual, 



= d. 
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and it follows that the Pauli equation in Lagrangian form 
becomes. 
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Henceforth, the derivatives are replaced by the covariant 
derivatives of a C/(l) x SU{2) gauge theory, where the 
SU{2) part takes care of the transport of spin. Surely, the 
second and especially the third term violate the SU{2) 
gauge invariance for the obvious reason that the non- 
Abelian 'gauge fields' A'^ are just proportional to the 

electromagnetic E and B fields. Notice that the second 
term just amounts to a small correction to the electro- 
magnetic part (third term) . The standard picture of how 
spins are precessing due to the spin-orbit coupling to ex- 
ternal electrical- and magnetic fields, pending the way 
they are moving through space can actually be taken as 
a literal cartoon of the parallel transport of non-Abelian 
charge in some fixed gauge potential! 

To be more precise, the SO problem does actually cor- 
respond with a particular gauge fix in the full SU{2) 
gauge theory. The electromagnetic fields have to obey 
the Maxwell equation. 
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and this in turn implies 
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Therefore, the SO problem is 'representative' for the 
SU(2) gauge theory in the Lorentz gauge and we do not 



have the choice of going to another gauge as the non- 
Abelian fields are expressed in terms of real electric and 
magnetic fields. This is a first new result. 

By varying the Lagrangian with respect to tp^ we ob- 
tain the Pauli equation in its standard Hamiltonian form. 
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where we leave the electromagnetic part implicit, antic- 
ipating that we will be interested to study the behavior 
of the quantum mechanical particles in fixed background 
electromagnetic field configurations. The wave function 
-0 can be written in the form. 



(10) 



with the probability density p, while 9 is the usual 
Abelian phase associated with the electromagnetic gauge 
fields. As imposed by the covariant derivatives, the 
SU{2) phase structure can be parametrised by the three 
non-Abelian phases ip'^, with the Pauli matrices r° act- 
ing on a reference spinor x- Hence, with regard to the 
wavefunction there is no difference whatever between the 
Pauli-problcm and genuine Yang-Mills quantum mechan- 
ics: this is all ruled by parallel transport. 

Let us now investigate in further detail how the Pauli 
equation transports spin-probability. This is in close con- 
tact with work in high-energy physics and we develop the 
theory along similar lines as Jackiwei al^. We intro- 
duce, however, a condensed matter inspired parametriza- 
tion that we perceive as instrumental towards laying bare 
the elegant meaning of the physics behind the equations. 

A key ingredient of our parametrization is the intro- 
duction of a non-Abelian phase velocity, an object occu- 
pying the adjoint together with the vector potentials. The 
equations in the remainder will involve time and space 
derivatives of 9, p and of the spin rotation operators 

e*"^"^"/^ . (11) 

Let us introduce the operator S°' as the non-Abelian 
charge at time t and at position r, as defined by the 
appropriate SU{2) rotation 
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(12) 



The temporal and spatial dependence arises through the 
non-Abelian phases (^"(i, r). The non-Abelian charges 
are, of course, SU{2) spin 1/2 operators: 



(13) 



It is illuminating to parametrize the derivatives of the 
spin rotation operators employing non-Abelian velocities 
u"- defined by, 



or 
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which are just the analogs of the usual Abelian phase 
velocity 



m m 



(15) 



These non-Abelian phase velocities represent the scale 
parameters for the propagation of spin probability in non- 
Abelian quantum mechanics, or either for the hydrody- 
namical flow of spin-superfluid. 

In addition we need the zeroth component of the ve- 
locity 



(16) 



being the time rate of change of the non-Abelian phase, 
amounting to a precise analog of the time derivative of 
the Abelian phase representing matter-density fluctua- 
tion, 



Mo = doe = -i—e^'^doe''' 



(17) 



It is straightforward to show that the definitions of the 
spin operators 5° , Eq. and the non-Abelian velocities 
Eq.'s fHl [TB)) . imply in combination. 
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It is easily checked that the definition of the phase 
velocity Eq. ([M]) implies the following identity. 



(19) 
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having as Abelian analogue, 

V X -u = . 



(20) 



as the latter controls vorticity, the former is in charge of 
the topology in the non-Abelian 'probability fluid'. It, 
however, acquires a truly quantum-hydrodynamical sta- 
tus in the rigid superfluid where it becomes an equation 
of algebraic topology. This equation is well known, both 
in gauge theory and in the theory of the '^He superfluids 
where it is known as the Mermin-Ho equatioiJ^. 



IV. SPIN TRANSPORT IN THE MESOSCOPIC 
REGIME 

Having defined the right variable, we can now go ahead 
with the quantum mechanics, finding transparent equa- 
tions for the non-Abelian probability transport. Given 
that this is about straight quantum mechanics, what fol- 
lows does bare relevance to coherent spin transport phe- 
nomena in the mesoscopic regime. We will actually de- 
rive some interesting results that reveal subtle caveats 
regarding mesoscopic spin transport. The punchline is 



that the Aharonov-Casher effect and related phenomena 
arc intrinsically fragile, requiring much more fine tun- 
ing in the experimental machinery than in the Abelian 
(Ahronov-Bohm) case. 

Recall the spinor definition Eq. fTOl) : together with the 
definitions of the phase velocity, it follows from the van- 
ishing of the imaginary part of the Pauli equation that, 



p[u A 

\ m 



(21) 



and this is nothing else than the non-Abelian continu- 
ity equation, imposing that probability is covariantly 
conserved. For non-Abelian parallel transport this is a 
weaker condition than for the simple Abelian case where 
the continuity equation implies a global conservation of 
mass, being in turn the condition for hydrodynamical 
degrees of freedom in the fluid context. Although lo- 
cally conserved, the non-Abelian charge is not globally 
conserved and this is the deep reason for the difficulties 
with associating a universal hydrodynamics to the non- 
Abelian fluids. The fluid dynamics will borrow this mo- 
tive directly from quantum mechanics where its meaning 
is straightforwardly isolated. 

Taking the trace over the non-Abelian labels in Eq. 
(|2ip results in the usual continuity equation for Abelian 
probability, in the spintronics context associated with the 
conservation of electrical charge. 
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where one recognizes the standard (Abelian) probability 
current. 
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From Abelian continuity and the full non-Abelian law Eq. 
pTjl it is directly seen that the non-Abelian velocities and 
vector potentials have to satisfy the following equations. 



n 



A" 



(24) 



and we recognize a divergence - the quantity inside the 
bracket is a conserved, current-like quantity. Notice that 
in this non-relativistic theory this equation contains only 
space like derivatives: it is a static constraint equation 
stating that the non-Abelian probability density should 
not change in time. The above is generally valid but it 
is instructive to now interpret this result in the Pauli- 
equation context. Using Eq.Q for the non Abelian vec- 
tor potentials, Eq ([M)) becomes. 



U^Ea 



(25) 



As a prelude to what is coming, we find that this actu- 
ally amounts to a statement about spin Hall probability 
currents. When the quantity on the r.h.s. would be zero, 
^aiiEi -I- V X A, the spin Hall equation 
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modulo an arbitrary curl and thus the spin Hall relation 
exhibits a "gauge invariance" . 

Let us complete this description of non-Abelian quan- 
tum mechanics by inspecting the real part of the Pauli 
equation in charge of the time evolution of the phase, 



do0 - eAa + u^S" - qA'^S" 



1 I m 

~h I 2" 

1 
2m 



-A + w°S'° 



1 A a 



2eg5M • A" + —A" • A'' 

4 



(26) 



h 

4m 



W (vp)2 



2p2 



Tracing out the non-Abelian sector we obtain the usual 
equation for the time rate of change of the Abelian phase, 
augmented by two SU{2) singlet terms on the r.h.s.. 



doe^eAo = 



4m 



2p2 



u A 

m 



2m 



T -A'' 



(27) 



Multiplying this equation by and tracing the non- 
Abelian labels we find. 



ug - gAg = fw - -A] ■ (u'' ~ ^A" 
n \ m / \ m 



(28) 



It is again instructive to consider the spin-orbit coupling 
interpretation, 



-ug = qBa - "T A, 

a \ m 



-^e,aiEi ) (29) 
m 



ignoring the spin orbit coupling this just amounts to Zee- 
man coupling. The second term on the right hand side 
is expressing that spin orbit coupling can generate uni- 
form magnetization, but this requires both matter cur- 
rent (first term) and a violation of the spin-Hall equation! 
As we have just seen such violations, if present, necessar- 
ily take the form of a curl. 

To appreciate further what these equations mean, let 
us consider an experiment of the Aharonov-Casher^ kind. 
The experiment consists of an electrical wire oriented, 
say, along the z-axis that is charged, and is therefore 
producing an electrical field Er in the radial direction in 
the xy plane. This wire is surrounded by a loop contain- 
ing mobile spin-carrying but electrically neutral particles 
(like neutrons or atoms ). Consider now the spins of the 
particles to be polarized along the z-direction and it is 
straightforward to demonstrate that the particles accu- 
mulate a holonomy ^ Er. It is easily seen that this cor- 
responds with a special case in the above formalism. By 
specializing to spins lying along the z-axis, only one com- 
ponent u^,Uq of the non-Abelian phase velocity u°',Uq 
has to be considered, and this reduces the problem to a 



U{1) parallel transport structure; this reduction is rather 
implicit in the standard treatment. 

Parametrise the current loop in terms of a radial (r) 
and azimuthal (0) direction. Insisting that the electrical 
field is entirely along r, while the spins are oriented along 
z and the current flows in the direction so that only 

7^ , Eq. dH]) reduces to (p(u^ - {q/m)Er)^ = 0. 
= pu^ corresponds with a spin probability current, 
and it follows that J| = {qp/m)Er + f{r,z) with / an 
arbitrary function of the vertical and radial coordinates: 
this is just the quantum-mechanical incarnation of the 
spin-Hall transport equation Eq. ([T])! For a very long 
wire in which all vertical coordinates are equivalent, the 
cylindrical symmetry imposes z independence, and since 
we are at fixed radius, / is a constant. In the case where 
the constant can dropped we have = d^9^ = {q/m)Er 
the phase accumulated by the particle by moving around 
the loop equals A9^ — § d<j)u^ — L{q/m)Er: this is just 
the Aharonov-Casher phase. There is the possibility that 
the Aharonov-Casher effect might not occur if physical 
conditions make the constant / nonzero. 

Inspecting the 'magnetization' equation Eq. (|29p . 
assuming there is no magnetic field while the parti- 
cle carries no electrical charge, = —{m/h)u- {u°- — 
{q/m)eiaiEi) = 0, given the conditions of the ideal 
Aharonov-Casher experiment. Henceforth, the spin cur- 
rents in the AC experiment do not give rise to magneti- 
zation. 

The standard AC effect appears to be an outcome of a 
rather special, in fact fine tuned experimental geometry, 
hiding the intricacies of the full non-Abelian situation 
expressed by our equations Eq. (|25l29p . As an example, 
let us consider the simple situation that, as before, the 
spins are polarized along the z-direction while the current 
flows along (j) such that only is non zero. However, we 
assume now a stray electrical field along the z-direction, 
and it follows from Eq. (|25|). 



d4> (p(i 



-Er) 



:uiE, 



(30) 



We thus see that if the field is not exactly radial, 
the nonradial parts will provide corrections to the spin 
Hall relation and more importantly will invalidate the 
Aharonov-Casher effect! This stray electrical field in 
the z-direction has an even simpler implication for the 
magnetization. Although no magnetization is induced in 
the z-direction, it follows from Eq. ([29)1 that this field 
will induce a magnetization in the radial direction since 
Uq = —u^{q/m)s^rzEz. This is finite since the matter 
phase current 7^ 0. 

From these simple examples it is clear that the non- 
Abelian nature of the mesoscopic spin transport under- 
lying the AC effect renders it to be a much less robust 
affair than its Abelian Aharonov Bohm counterpart. In 
the standard treatment these subtleties are worked under 
the rug and it would be quite worthwhile to revisit this 
physics in detail, both experimentally and theoretically, 
to find out if there are further surprises. This is however 
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not the aim of this paper. The general message is that 
even in this rather weU behaved mesoscopic regime al- 
ready finds the first signs of the fragility of non-Abelian 
transport. On the one hand, this will turn out to be- 
come lethal in the classical regime, while on the other 
hand we will demonstrate that the coherent transport 
structures highlighted in this section will acquire hydro- 
dynamical robustness when combined with the rigidity 
of non-Abelian superfluid order. 



derivatives for the normal derivatives 

r = 

= JS" 



ih 
2m 



p 

4 \ m 

fa I ja 



A" 



(35) 



(36) 



where the gauged version of the non-coherent and coher- 
ent currents are respectively, 



V. SPIN CURRENTS ARE ONLY 
COVARIANTLY CONSERVED. 

It might seem odd that the quantum equations of the 
previous section did not have any resemblance to a con- 
tinuity equation associated with the conservation of spin 
density. To make further progress in our pursuit to de- 
scribe macroscopic spin hydrodynamics an equation of 
this kind is required, and it is actually straightforward 
to derive using a different strategy (see also Jackiw et 

Let us define a spin density operator, 



and a spin current operator, 



J 



ih 
3nc 



V'^YV^-(VV')^yV 



Jc 



(31) 



(32) 



We observe that the spin current operator can be written 
as a sum of two contributions. The first piece can be 
written as 



(33) 



It factors in the phase velocity associated with 
the Abelian mass current u times the non-Abelian 
charge/spin density E° carried around by the mass cur- 
rent. This 'non-coherent' (relative to spin) current is 
according to the simple classical intuition of what a spin 
current is: particles flow with a velocity u and every par- 
ticle carries around a spin. The less intuitive, 'coherent' 
contribution to the spin current needs entanglement of 
the spins, 



(34) 



and this is just the current associated with the non- 
Abelian phase velocity u"" already highlighted in the pre- 
vious section. 

The above expressions for the non-Abelian currents are 
of relevance to the 'neutral' spin fluids, but we have to 
deal with the gauged currents, for instance because of 
SO-coupling. Obviously we have to substitute covariant 



ja J Qa 

•Jnc — ■J ^ 



J n — ~ 



^1" 



(37) 
(38) 



with the Abelian (mass) current J given by Eq. 

It is a textbook exercise to demonstrate that the follow- 
ing 'continuity' equations holds for a Hamiltonian char- 
acterized by covariant derivatives (like the Pauli Hamil- 
tonian), 



(39) 



with the usual non-Abelian covariant derivatives of 
vector-fields. 



(40) 



Eq. (|39p has the structure of a continuity equation, ex- 
cept that the derivatives are replaced by covariant deriva- 
tives. It is well known^ that in the non-Abelian case 
such covariant 'conservation' laws fall short of being real 
conservation laws of the kind encountered in the Abelian 
theory. Although they impose a local continuity, they 
fail with regard to global conservation because they do 
not correspond with total derivatives. This is easily seen 
by rewriting Eq. p9)) as 



9 abc^ . J'' 



(41) 



The above is standard lore. However, using the result 
Eq. ([M]) from the previous section, we can obtain a bit 
more insight in the special nature of the phase coher- 
ent spin current, Eq. Eq. ([M]) can be written in 
covariant form as 



involving only the space components and therefore 



(42) 



(43) 



Since E'' is spin density, it follows rather surprisingly 
that the coherent part of the spin current cannot give 
rise to spin accumulation] Spin accumulation is entirely 
due to the non-coherent part of the current. Anticipating 
what is coming, the currents in the spin superfluid are 
entirely of the coherent type and this 'non-accumulation 
theorem' stresses the rather elusive character of these 
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spin supercurrents: they are so 'unmagnetic' in character 
that they are even not capable of causing magnetization 
when they come to a standstill due to the presence of a 
barrier! 

As a caveat, from the definitions of the coherent- and 
non-coherent spin currents the following equations can 
be derived 

p (V X J^c) - 4^e"''Vl, X + x 

(44) 

- Ipe'^"^!" ■ ■ (45) 

From these equations it follows that the coherent currents 
actually do influence the way that the incoherent currents 
do accumulate magnetization, but only indirectly. Simi- 
larly, using the divergence of the Abclian covariant spin 
current together with the covariant conservation law, we 
obtain the time rate of precession of the local spin density 

doT," = 1^5" + 4^e°^'^ • Jnc " le^'^^oS" • (46) 

demonstrating that this is influenced by the presence of 
coherent- and incoherent currents flowing in orthogonal 
non-Abelian directions. 

This equation forms the starting point of the discus- 
sion of the (lack of) hydrodynamics of the classical non- 
Abelian/spin fluid. 

VI. PARTICLE BASED NON-ABELIAN 
HYDRODYNAMICS, OR THE CLASSICAL SPIN 
FLUID. 

We have now arrived at a point that we can start to 
address the core-business of this paper: what can be 
said about the collective flow properties of large assem- 
blies of interacting particles carrying spin or either non- 
Abelian charge? In other words, what is the meaning of 
spin- or non-Abelian hydrodynamics? The answer is: if 
there is no order-parameter protecting the non-Abelian 
phase coherence on macroscopic scales spin flow is non- 
hydrodynamical, i.e. macroscopic flow of spins does not 
even exist. 

The absence of order parameter rigidity means that we 
are considering classical spin fluids as they are realized 
at higher temperatures, i.e. away from the mesoscopic 
regime of the previous section and the superfluids ad- 
dressed in Section IIXI The lack of hydrodynamics is well 
understood in the spintronics community: after generat- 
ing a spin current is just disappears after a time called 
the spin-relaxation time. This time depends of the ef- 
fective spin-orbit coupling strength in the material but 
it will not exceed in even the most favorable cases the 
nanosecond regime, or the micron length scale. Surely, 
this is a major (if not fundamental) obstacle for the use of 



spin currents for electronic switching purposes. Although 
spin currents are intrinsically less dissipative than elec- 
trical currents it takes a lot of energy to replenish these 
currents, rendering spintronic circuitry as rather useless 
as competitors for Intel chips. 

Although this problem seems not to be widely known 
in corporate head quarters, or either government fund- 
ing agencies, it is well understood in the scientific com- 
munity. This seems to be a different story in the 
community devoted to the understanding of the quark- 
gluon plasma's produced at the heavy ion collider at 
Brookhaven. In these collisions a 'non-Abelian fire ball' 
is generated, governed by high temperature quark-gluon 
dynamics: the temperatures reached in these fireballs ex- 
ceed the confinement scale. To understand what is hap- 
pening one of course needs a hydrodynamical description 
where especially the fate of color (non-Abelian) currents 
is important. It seems that the theoretical mainstream 
in this pursuit is preoccupied by constructing Boltzmann 
type transport equations. Remarkably, it does not seem 
to be widely understood that one first needs a hydro- 
dynamical description, before one can attempt to calcu- 
late the numbers governing the hydrodynamics from mi- 
croscopic principle by employing kinetic equations (quite 
questionable by itself given the strongly interacting na- 
ture of the quark-gluon plasma). The description of the 
color currents in the quark-gluon plasma is suffering from 
a fatal flaw: because of the lack of a hydrodynamical con- 
servation law there is no hydrodynamical description of 
color transport. 

The above statements are not at all original in this re- 
gard: this case is forcefully made in the work by Jackiw 
and coworkersi^iii dealing with non-Abelian 'hydrody- 
namics'. It might be less obvious, however, that pre- 
cisely the same physical principles are at work in the 
spin-currents of spintronics: spintronics can be viewed 
in this regard as 'analogous system' for the study of the 
dynamics of quark-gluon plasma's. The reason for the 
analogy to be precise is that the reasons for the failure of 
hydrodynamics reside in the parallel transport structure 
of the matter fields, and the fact that the 'gauge fields' 
of spintronics are in 'fixed frame' is irrelevant for this 
particular issue. 

The discussion by Jackiw et al. of classical ('particle 
based') non-Abelian 'hydrodynamics' starts with the co- 
variant conservation law we re-derived in the previous 
section, Eq. This is still a microscopic equation 

describing the quantum physics of a single particle and 
a coarse graining procedure has to be specified in order 
to arrive at a macroscopic continuity equation. Rest- 
ing on the knowledge about the Abelian case this coarse 
graining procedure is unambiguous when we are inter- 
ested in the (effective) high temperature limit. The nov- 
elty as compared the Abelian case is the existence of the 
coherent current Jq expressing the transport of the en- 
tanglement associated with non-Abelian character of the 
charge; Abelian theory is special in this regard because 
there is no room for this kind of entanglement. By defini- 
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tion, in the classical limit quantum entanglement cannot 
be transported over macroscopic distances and this im- 
plies that the expectation value (J^) cannot enter the 
macroscopic fluid equations. Although not stated explic- 
itly by Jackiw et al, this particular physical assumption 
(or definition) is the crucial piece for what follows - the 
coherent current will acquire (quantum) hydrodynamic 
status when protected by the order parameter in the spin- 
superfluids. 

What remains is the non-coherent part, governed by 
the pseudo-continuity equation Eq. (j43p . Let us first con- 
sider the case that the non-Abelian fields are absent (e.g., 
no spin-orbit coupling) and the hydrodynamical status of 
the equation is immediately obvious through the Ehrcn- 
fest theorem. The quantity T,"- (p<S'°) becomes just 
the macroscopic magnetization (or non-Abelian charge 
density) that can be written as nQ, i.e. the macro- 
scopic particle density n — (p) times their average spin 
Q = (S). Similarly, the Abelian phase current pu turns 
into the hydrodynamical current nv where v is the veloc- 
ity associated with the macroscopic 'element of fluid'. In 
terms of these macroscopic quantities, the l.h.s. of Eq. 
([29|l just expresses the hydrodynamical conservation of 
uniform magnetization in the absence of spin-orbit cou- 
pling. In the presence of spin orbit coupling (or gluons) 
the r.h.s. is no longer zero and, henceforth, uniform mag- 
netization/color charge is no longer conserved. 

Upon inserting these expectation values in Eq.'s (P^ . 
P5)) one obtains the equations governing classical non- 
Abelian fluid flow, 



Translating the gauge fields into the physical electromag- 
netic fields of the Pauli equation, Eq. ([^5]) becomes. 



dtn + V • (nv) — 



(47) 
(48) 



Eq. (|47p expresses the usual continuity equation associ- 
ated with (Abelian) mass density. Eq. ([48)1 is the novelty, 
reflecting the non-Abelian parallel transport structure, 
rendering the substantial time derivative of the magne- 
tiziation/color charge to become dependent on the color 
charge itself in the presence of the non-Abelian gauge 
fields. To obtain a full set of hydrodynamical equations, 
one needs in addition a 'force' (Navier-Stokes) equation 
expressing how the Abelian current nv accelerates in the 
presence of external forces, viscosity, etcetera. For our 
present purposes, this is of secondary interest and we re- 
fer to Jackiw et aL^-'^^for its form in the case of a perfect 
(Euler) Yang-Mills fluid. 

Jackiw e< al. coined the name 'Fluid- Wong Equa- 
tions' for this set of equations governing classical non- 
Abelian fluid flow. These would describe a hydrody- 
namics that would be qualitatively similar to the usual 
Abelian magneto-hydrodynamics associated with electro- 
magnetic plasma's were it not for Eq. this expres- 
sion shows that the color charge becomes itself depen- 
dent on the flow. This unpleasant fact renders the non- 
Abelian flow to become non-hydrodynamical. 

We perceive it as quite instructive to consider what 
this means in the spintronics interpretation of the above. 



dtQ^' + v-VQ'^ - ( 



cB + v X E 



X Q 



(49) 



where Q{r) has now the interpretation of the uniform 
magnetization associated with the fluid element at po- 
sition r. The flrst term on the r.h.s. is just expressing 
that the magnetization will have a precession rate in the 
comoving frame, proportional to the external magnetic 
field B. However, in the presence of spin-orbit coupling 
(second term) this rate will also become dependent on 
the velocity of the fluid element itself when an electrical 
fleld E is present with a component at a right angle both 
to the direction of the velocity v and the magnetization 
itself. This velocity dependence wrecks the hydrodynam- 
ics. 

The standard treatments in terms of Boltzmann equa- 
tions lay much emphasis on quenched disorder, destroy- 
ing momentum conservation. To an extent this is obscur- 
ing the real issues, and let us instead focus on the truly 
hydrodynamical flows associated with the Galilean con- 
tinuum. For a given hydrodynamical flow pattern, elec- 
tromagnetic field configuration and initial configuration 
of the magnetization, Eq. (|49)) determines the evolution 
of the magnetization. Let us consider two elementary ex- 
amples. In both cases we consider a Rashba-lik^i elec- 
tromagnetic field configuration: consider fiow patterns in 
the xy directions and a uniform electrical fleld along the 
z direction while B — 0. 

a. Laminar flow 

Consider a smooth, non-turbulent laminar flow pattern 
in a 'spin-fluid tube' realized under the condition that the 
Reynold's number associated with the mass flow is small. 
Imagine that the fluid elements entering the tube on the 
far left have their magnetization Q oriented in the same 
direction (Fig. [2]). Assume first that the velocity v is 
uniform inside the tube and it follows directly from Eq. 

that the Q's will precess with a uniform rate when 
the fluid elements move trough the tube. Assuming that 
the fluid elements arriving at the entry of the tube have 
the same orientation at all times, the result is that an ob- 
server in the lab frame will measure a static 'spin spiral' 
in the tube, see Fig. [31 At flrst sight this looks like the 
spiral spin structures responsible for the ferroelectricity 
in the multiferroics but this is actually misleading: as we 
will see in Section VII these are actually associated with 
localized particles (i.e. no Abelian flow) while they are 
rooted instead in the entanglement current. We leave it 
as an excercise for the reader to demonstrate that the 
spiral pattern actually will not change when the flow in 
the tube acquires a typical laminar, non-uniform velocity 
distribution, with the velocities vanishing at the walls. 

b. Turbulent flow 

Let us now consider the case that the fluid is moving 
much faster, such that downstream of an obstruction in 
the flow turbulence arises in the matter current. In Fig- 
ure m we have indicated a typical stream line showing 
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Non-Abelian charge 



laminar flow 




FIG. 2: Laminar flow of a classical spin fluid in an electric 
field. The fluid elements (blue) carry non-Abelian charge, the 
red arrows indicating the spin direction. The flow lines are 
directed to the right, and the electric field is pointing outwards 
of the paper. Due to Eq. (|49|l . the spin precesses as indicated. 

t / ' \ \ / ^ 

\ / ■ \ i / ■ 

' \ i / 

\ / ' \ 1 / 

FIG. 3: The laminar flow of a parallel transported spin cur- 
rent. Figure [2j can also be viewed as a static spin spiral mag- 
net. 



with regard to the possibility that an elegant description 
will be found, in a way mirroring the state of spintronics. 
We will instead continue now with our exposition of the 
remarkable fact that the rigidity associated with order 
parameters is not only simplifying the hydrodynamics 
(as in the Abelian case) but even making it possible for 
hydrodynamics to exist! 




FIG. 4: Turbulent spin flow around an obstruction in an 
electric field. It is seen that only the "mass" is conserved. 
The change in spin direction after one precession around the 
obstruction causes a spin sink. Hence it is precisely the par- 
allel transport, or the covariant conservation, which destroys 
hydrodynamic conservation for non-Abelian charge. 



that the flow is now characterized by a finite vorticity 
in the region behind the obstruction. Let us now repeat 
the exercise, assuming that fluid elements arrive at the 
obstruction with aligned magnetization vectors. Follow- 
ing a fluid element when it traverses the region with finite 
circulation it is immediately obvious that even for a fixed 
precession rate the non-Abelian charge/magnetization be- 
comes multivalued when it has travelled around the vor- 
tex! Henceforth, at long times the magnetization will 
average away and the spin current actually disappears 
at the 'sink' associated with the rotational Abelian flow. 
This elementary example highlights the essence of the 
problem dealing with non-Abelian 'hydrodynamics': the 
covariant conservation principle underlying everything 
is good enough to ensure a local conservation of non- 
Abelian charge so that one can reliably predict how the 
spin current evolves over infinitesimal times and dis- 
tances. However, it fails to impose a global conservation. 
This is neatly illustrated in this simple hydrodynamical 
example: at the moment the mass fiow becomes topo- 
logically non-trivial it is no longer possible to construct 
globally consistent non-Abelian flow patterns with the 
consequence that the spin currents just disappear. 

Although obscured by irrelevant details, the above mo- 
tive has been recognized in the literature on spin flow 
in semiconductors where it is known as D'yakonov-Perel 
spin relaxation^^, responsible for the longitudinal (Ti) 
spin relaxation time. We hope that the analogy with 
spin-transport in solids is helpful for the community that 
is trying to flnd out what is actually going on in the 
quark-gluon flreballs. Because one has to deal eventually 
with the absence of hydrodynamics we are pessimistic 



VII. ELECTRODYNAMICS OF SPIN-ORBIT 
COUPLED SYSTEMS 

Before we address the interesting and novel effects in 
multiferroics and spin superfluids, we pause to obtain the 
electrodynamics of spin orbit coupled systems. From the 
Pauli Maxwell Lagrangian ^ we see that the spin cur- 
rent couples directly to the electric fleld and will thus 
act as a source for electric fields. In order to see how 
this comes about let us obtain the electrodynamics of 
a spin-orbit coupled system. We presuppose the usual 
definition of electromagnetic fields in terms of gauge po- 
tentials, which implies the Maxwell equations 

V ■ B = 0, V X E + doB = 0. (50) 

If we vary the Lagrangian with respect to the scalar elec- 
tromagnetic potential, we obtain 

d,E,^ATrqe,ai{x^dai'x) (51) 

where we suppose that the charge sources are cancelled 
by the background ionic lattice of the material or that we 
have a neutral system. This term is extremely interesting 
because it says that the "curl" of spin currents are sources 
for electric fields. In fact, the electric field equation is 
nothing but the usual Maxwell equation for the electric 
displacement V • 13 = where D = E + AnP with 

P^ - -e^aiX^Jl'x ■ (52) 

The spin current acts as an electrical polarization for the 
material. The physical origin of this polarization is rela- 
tivistic. In the local frame the moving spins in the current 
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produce a magnetic field as they are magnetic moments. 
After a Lorentz transformation to the lab frame, part of 
this field becomes electric. On the other hand, it can be 
shown that V-P = unless the spin current has singular- 
ities. Thus, in the absence of singularities spin currents 
cannot create electric fields. 

Varying the Lagrangian ^ with respect to the vector 
potential we obtain 

(\/xBy^ An Jem - 47r X + OqE^ 

- A7Tqela^^o (x^Tx) (53) 
= An Jem ~An(v X gE^ + doD, . 

The first term on the right hand side contains the usual 
electromagnetic current 

Jem = Anep{u,+u'}x^S''x) (54) 

which includes the motion of particles due to the ad- 
vance of the Abelian and the non-Abelian phases. The 
term containing the non-Abelian velocity (the coherent 
spin current) in this electromagnetic current will only 
contribute when there is magnetic order (5*°) ^ 0. The 
second term is conventional since it is the curl of the mag- 
netization which generates magnetic fields. The third is 
the Maxwell displacement current in accordance with our 
identification of the electrical polarization caused by the 
spin current. 

VIII. SPIN HYDRODYNAMICS RISING FROM 
THE ASHES I: THE SPIRAL MAGNETS. 

Recently the research in multiferroics has revived. This 
refers to materials that are at the same time ferroelectric 
and ferromagnetic, while both order parameters are cou- 
pled. The physics underlying this phenomenon goes back 
to the days of Lifshitz and Landau^^. Just from consid- 
erations regarding the allowed invariants in the free en- 
ergy it is straightforward to find out that when a crystals 
lacks an inversion center (i.e., there is a net internal elec- 
tric field) spin-spin interactions should exist giving rise 
to a spiral modulation of the spins (helicoidal magnets). 
The modern twist of this argument is2.: the spin spiral 
can be caused by magnetic frustration as well, and it now 
acts as a cause (instead of effect) for an induced ferro- 
electric polarization. Regarding the microscopic origin of 
these effects, two mechanisms have been identified. The 
first one is called 'exchange striction' and is based on the 
idea that spin-phonon interactions of the kind familiar 
from spin-Peierls physics give rise to a deformation of the 
crystal structure when the spin-spiral order is present, 
and these can break inversion symmetry^. The second 
mechanism is of direct relevance to the present subject 
matter. As we already explained in the previous section, 
a spiral in the spin-density can be viewed at the same 
time as a spin current. In the presence of the magnetic 



order parameter this spin current acquires rigidity (like 
a supercurrent) and therefore it can impose its will on 
the 'gauge' fields. In the spin-orbital coupling case, the 
'gauge' field of relevance is the physical electrical field, 
and henceforth the 'automatic' spin currents associated 
with the spiral magnet induce an electrical field via the 
spin-orbit coupling, rendering the substance to become a 
ferroelectric'^. 

This substance matter is rather well understood^ and 
the primary aim of this section is to explain how these 
'spiral magnet' spin currents fit into the greater picture 
of spin-hydrodynamics in general. Viewed from this gen- 
eral perspective they are quite interesting: they belong 
to a category of non-Abelian hydrodynamical phenom- 
ena having no analogy in the Abelian universe. On the 
one hand these currents are spontaneous and truly non- 
dissipative and in this regard they are like Abelian super- 
currents. They should not be confused with the Froehlich 
'super' currents associated with (Abehan) charge density 
waves: these require a time dependence of the density or- 
der parameter (i.e., the density wave is sliding) while the 
spiral magnet currents flow also when the non-Abelian 
density (the spiral) is static. This belies their origin in 
the coherent non-Abelian phase current just as in the 
spin-superfluids, or either the non-Abelian Higgs phase. 

An important property of the static coherent spin cur- 
rents of the spin spirals is that vortex textures in the 
spin background become sources of electrical charge in 
the presence of spin-orbit coupling, as first observed by 
Mostovoy^ . Anticipating the discussion of the SO cou- 
pled spin superfluid in the next sections, a major differ- 
ence between those and the multiferroics is that in the 
former the phase coherent spin fluid can quantize the elec- 
trical line charge but not cause electrical charge because 
of the important difference that such a current can not 
originate spontaneously in the spin superfluid because 
it needs to be created by an electric field. It can trap 
charge because being a supercurrent it does not decay if 
the battery that creates the electric field is removed. 

Last but not least, the spiral magnet currents offer a 
minimal context to illustrate the most fundamental fea- 
ture of non-Abelian hydrodynamics: the rigidity of the 
order parameter is capable of restoring hydrodynamical 
degrees of freedom that are absent in the 'normal' fluid 
at high temperature. This is so simple that we can ex- 
plain it in one sentence. One directly recognizes the XY 
spin vortex in the turbulent flow of Fig. [4l but in the 
presence of spin density order the 'spiral' spin pattern 
associated with the vortex has to be single valued, and 
this in turns renders the spin current to be single valued: 
spin currents do not get lost in the ordered magnet! 

To become more explicit, let us rederive Mostovoy's 
result in the language of this paper, by considering an 
ordered XY-magnet with an order parameter that is the 
expectation value of the local spin operator 

{S^ + iSy) = Se'^ (55) 
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In general a spin state of an XF-niagnet is given by 

n 3(^)1 T> (56) 

lattice sites 

where we specialize to spin 1/2 for explicitness, but sim- 
ilar results hold for larger spin. | t) is a spinor in the 
+z direction and g(x) is an SU (2) rotation matrix in the 
xy-plane: 

g{x) = e*«(^)^-/2 (57) 

where is the Pauli matrix in the z-direction. The or- 
dered ground state of the uniform XY-magnet requires 
that 0{x) and hence g{x) are independent of x. Besides 
the ground state, XF-magnets have excited metastable 
states corresponding to XY spin vortices. These are eas- 
ily constructed by choosing 



ncf) , n integer , 



arctan 



(f) 



(58) 



Now we can compute the spin current in this state. The 
coherent spin current is given by 



Jc = i^u'^ 



. tip 
2m 



(59) 



For our case 



1 



9 = 2 ^'^^cosw 

9 y5 = 2 



(60) 



we have the appropriate 0(2) or ?7(1) rotation. We also 
have for the vortex 9 — mp 



%Lpr^ 12 f a 



{t ,t |( 



nhp 
Am 



(61) 



(V(^)(5''^ 



According to the results in the previous section, spin cur- 
rents alter the electrodynamics via Gauss' law, 



(62) 



where q measures the coupling between spin currents and 
electric fields via spin orbit coupling. Hence, using that 
for = arctan(y/x). 



V X V0 2nS^^\f) 
we find for the spin current of the vortex, 

d,E,^2^nq^6^^Hr). 



(63) 



(64) 



Therefore spin vortices in XF-magnets produce electric 
fields! 



IX. SPIN HYDRODYNAMICS RISING FROM 
THE ASHES II: THE SPIN SUPERFLUIDS 

Even without knowing a proper physical example of a 
spin-orbit coupled spin-superfluid one can construct its 
order parameter theory using the general principles dis- 
covered by Ginzburg and Landau. One imagines a con- 
densate formed from electrically neutron bosons carrying 
SU{2) spin triplet quantum numbers. This condensate 
is characterized by a spinorial order parameter. 



^(i9+i(p"r"/2) 



X 



(65) 



where j^*! is the order parameter amplitude, nonzero in 
the superfluid state, while 9 is the usual U{1) phase asso- 
ciated with number, while the three non-Abelian phases 
with the Pauli matrices t° acting on a reference 
spinor x keep track of the SU{2) phase structure. Ac- 
cording to the Ginzburg-Landau recipe, the free energy 
of the system should be composed of scalars constructed 
from while the gradient structure should be of the same 
covariant form as for the microscopic problem - parallel 
transport is marginal under renormalization. Henceforth, 
we can directly write down the Ginzburg-Landau free en- 
ergy density for the spin superfluid in the presence of spin 
orbit coupling. 



fa ^ 
2m 



1 

8^' 



2m ^ 4 ^ 



(66) 



We now specialize to the deeply non-relativistic case 
where the time derivatives can be ignored, while we con- 
sider electrically neutral particles (e = 0) so that the EM 
gauge fields drop out from the covariant derivatives. 

Well below the superfluid transition the amplitude l^fj 
is finite and frozen and one can construct a London-type 
action. Using the formulas in the appendix we obtain 
that 



r ■ 

*^sr)i 



spin-vcl 



m 



8m 



(67) 

Using the spin identities defined in Section |Vl this can 
be rewritten as 



-'spin-vcl 



c\ ja ja £) ja ja 



.J—A^'-A''. 
8m 



(68) 

We see that the Ginzburg-Landau action is a sum of the 
spin coherent and non-coherent squared currents. The 
spin noncoherent part has to do with mass or U{1) cur- 
rents, but since the particles carry spin they provide a 
spin current only if (5*°) ^ 0, requiring a net magne- 
tization. The coherent part is a bona fide spin current 
originating in the coherent advance of the non-Abelian 
phase associated with the spin direction. 
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In order to make contact with the Hehum literatureiS- 
we will write our spin operators and the coherent spin 
currents in terms of SO{3) rotation matrices via 



R%iv) Y - e-'^^"^"/2 ^ 6^'^"-°/^ (69) 



with R°fj{ip) an 5*0(3) rotation matrix around the vector 
ip by an angle |</?|, we obtain that the spin operator is a 
local 50(3) rotation of the Pauli matrices 



(70) 



In terms of the rotation operators, the spin velocities 
related to advance of the non-Abelian phase are 



TO 



It is also easily seen that 



(71) 



(72) 



If we look at the expressions for it" and in terms of 
the spin rotation matrix for the spin-orbit coupled spin 
superfluid, Eq.'s ([TTl [72]), we recognize these to be the 
exact analogues of the spin velocity and spin angular ve- 
locity of ^He-B (jllip reproduced in Section IXI Al We 



define g through 
that is 



-1 ' a 



9 ^y.? = Sa 



g = e 



(73) 



(74) 



which is an SU(2) group element. We now have the spin 
velocities and angular velocities expressed as 



= -iT:T{Sag ^dig] = -iTY^g ^y^^s} 
^ -iTr [Sag^^dog] = -iTr |.g"^y9o5| 



(75) 



The first is proportional to the coherent spin current and 
the second to the effective magnetization. If we define 
the spin superfiuid density via 



1 2 
P = -^XBC 

T 



(76) 



we have the following Lagrangian that describes the low 
energy spin physics, written in a way that is quite anal- 
ogous to that of ^He-Bi^, 



1 

2^2 



„2 2 



he 



— (E^ - B^) 



(77) 



From this Lagrangian we obtain the spin equations of 
motion for the spin superfiuid by varying with respect to 
the non-Abelian phase 



We evaluate 



^0 



" dL ■ 




■ dL ' 


[d{dog)\ 




[d{d,g)\ 



dL 
dg 
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(78) 



dL dg ^ duia dL dg ^ duoai dL 



dg dg dg~^ duo a dg dg~^ duoo 
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dL 
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dL 
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(79) 
(80) 

(81) 
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which yields the rather formidable equation of motion 
= do 

H-a- 



■ -1 '''a 2 I 7^ 



2^1 



(82) 



*5 



After some straightforward algebra this equation reduces 
to the fairly simple equation 



do {uJa + iBa) ~ C^di UJa 



2^ 

he 



f-aikEk 1 — 



(83) 



The solution of this equation of motion gives the SU(2) 
group element 5 as a function of space and time, and the 
spin velocities and angular velocities can be determined. 

Similarly, by varying the Lagrangian (|83|) with respect 
to the electromagnetic potentials, we obtain the Maxwell 
equations for the electromagnetic fields "created" by the 
spin velocities and angular velocities. 



^ /2c^ 



(V X = -47r X ^XB^^ 

2c/i 



(84) 



(85) 



We like to draw the reader's attention to the fact that 
Mineev and Volovik derived these results already in the 
seventies^ in the context of '^He-B. We show here that 
these hold in the general case of an SU{2) spin super- 
fluid, and will demonstrate in section IXI CI that similar 
equations can be derived for the case of superfluid '^He-A 
as well. 



X. CHARGE TRAPPING BY SPIN 
SUPERFLUIDS 

We now go back to the trick of charge trapping in su- 
perfluids we used previously to wet your appetite. How 
docs this magic trick work? At the heart of our idea lies 
the spin vortex solution. Let us first briefly sketch the 
argument, and then prove it. The straight wire causes 
an electric field of 



E = 



X 



2neor 



(86) 



where f is a radial unit vector in the xy plane perpendic- 
ular to the cylinder axis z. The azimuthal angle is if. We 
now need to determine the electric field in the superfiuid 



region. Because of the symmetry of the problem, this 
electric field will be radial. Lets call it E^. This electric 
field will drive a spin current, which will be a source of 
electric field itself if it has a singularity that will lie on 
the wire because of the radial symmetry. The symmetry 
of the problem suggests that the spins will be polarized 
along the axis of the cylinder. By solving the equations 
of motion in the presence of an electric field and no mag- 
netic field, we obtain that when the spin current and spin 
angular velocity satisfy the Spin Hall relation for spin di- 
rection a — z 



(87) 



with the magnetic moment of the He-atoms 

me 

M = g MB, 

rriHe 

whereas the other spin superfluid velocities vanish. Since 
the electric fields do not depend on the z-coordinate and 
only have a radial component, the equations of motion 
Eq. (|83p are satisfied. In our case, written in cylindrical 
coordinates. 



2^ 



(89) 



We see that the electric field leads to a spin vortex, i.e., 
z-polarised spins flowing around the wire. This is noth- 
ing different from vortices in Bose superfiuids induced by 
rotation. This might cause some concern as we have an 
SU{2) superfluid while vortices are topological defects 
associated with U{1). Why is this spin vortex topologi- 
cally stable? This has everything to do with the fact that 
we are not dealing with a real gauge theory but that our 
'gauge' fields are in fact physical. In a literal sense, the 
topology is 'hard wired' by the fact that we have put 
the wire inside the cylinder: the electrical field is forced 
by the experimentalist to carry a vortex topology, and 
this topology is via the covariant derivatives imposed on 
the spin current - were it a real (unphysical) gauge field, 
it has to sort this out by itself and the outcome would 
be the usual 't Hooft-Polyakov monopole. There is a 
neat mathematical way of saying the same thing. Gauge 
theories coupled to matter are known to mathematicians 
as bundle theories. One way to classify them is by us- 
ing Chern classe a^^i^^ . The Chern classes do not depend 
on the gauge chosen, or the configuration of the matter 
fields, but are a property of the bundle. The ramification 
is that if the topology of the gauge field is cylindrical, the 
matter field has cylindrical topology as well. 

The stability of the vortex can also be checked by 
demonstrating that a vortex centered on the wire, with a 
spin direction parallel to this wire, does satisfy the equa- 
tions of motion we derived in Section IX, while such a 
solution is an energy minimum. From the Lagrangian in 
the previous section it follows that the momentum con- 
jugate to the non-Abelian phase is 



n = 



XBC. 

272 



^UJe,^eo.^kEk] + ^E^ . (90) 
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When the vortex solution ,and thereby the Spin Hall re- 
lation is valid, the energy density becomes, 



If there is no vortex we have energy density 



1 



'^no-vortcx — ^ ^ 

Stt 



(91) 



(92) 



which is bigger than the energy density TLsh correspond- 
ing to a vortex present and thus the solution with the vor- 
tex is favored. If we have a vortex solution and perturb 
around by SijJai the energy changes by 



51-L 



(93) 



which is a positive quantity and we see that the vortex 
solution is stable against perturbations as they increase 
the energy of the system. We can rephrase the above 
reasoning in a more sophisticated way: the cylindrical 
topology of the fixed-frame gauge fields imposes the same 
vortex-type topology on the matter field, because of the 
parallel transport structure originating from spin-orbit 
coupling! 

The vortex topology can be classified by winding num- 
bers. Indeed, from the definition of the spin supercurrent 
in chapter ITXl we have 



to. 



-ve. 



(94) 



Therefore the spin current must satisfy the quantization 
condition 



• dl = 27riV 



(95) 



when we integrate around the cylinder where N is an 
integer. This quantisation is not quite shocking, since 
any order parameter theory has this condition. However, 
bearing in mind the magnetic flux trapping in supercon- 
ductors, it is interesting to integrate the spin current af- 
ter substituting the spin-Hall equation. By Gauss' law, 
one obtains that the very same phase velocity integral 
becomes 



LOz ■ dl2Tr fJ-o^- 

rriHe 



(96) 



In other words, the charge density is quantised in units 
of 

A = iVAo = = 2.6 X lO-^C/m!. (97) 

in the specific case of ^HeThis is of course a very large 
line-charge density, and this is of course rooted in the fact 
that this quantum is 'dual' to the tiny spin orbit coupling 
of helium, in the same way that the flux quantum in 
superconductors is inversely proportional to the electrical 



charge. In he next section we will show that this huge 
required electrical charge is detremental to any attempt 
to realize such an experiment employing a substance like 
helium. 

This experiment is the rigid realisation of the 
Aharonov-Casher phase^, for which our application is 
inspired by Balatskii and Altshuler~. The rigidity is 
provided by the superfluid density, forcing the winding 
number to be integer. Our idea is actually the spin super- 
fluid analogue of the flux trapping with superconducting 
rings. The quantization of magnetic flux is provided by 
the screening of electromagnetic fields, causing vanishing 
total superconducting current. The latter, being defined 
covariantly, consists of a U{1) superfluid velocity and a 
gauge fleld. Calculating the line integral 



= (/) Jrdx, ^ id: 



Ajdxi — 27m — 



(98) 



leading to the flux quantisation condition. In the above 
argument, the gauge fields Ai have dynamics, leading to 
screening of the Ai in the superconducting ring. 

In our case, the gauge fields are fixed by the electro- 
magnetic fields, such that there cannot be screening ef- 
fects. Still, the spin-Hall equations, which solve the equa- 
tions of motion (j83p , lead to a vanishing superconducting 
current. The gauge fields, being unscreened, play now a 
quite different role: these are necessary to force the topol- 
ogy of the superfluid order parameter to be U{1). The 
result is the same: quantisation of electric flux, deter- 
mined by the charge on the wire. 

Charge trapping in spin superfluids and in magnets 
both originate from the coherent part of the spin current. 
In this sense, there is not too much difference between 
the two effects. On the other hand, there is a subtle, 
but important distinction. For magnets there is no need 
for electric fields to impose the supercurrent, since they 
are wired in by the magnetic order. In contrast in the 
spin superfluids, an electric fleld is necessary to create a 
coherent spin current since there is no magnetisation, and 
in this sense the spin superfluids cannot create electrical 
charge, while magnets can. 

The question which surely is nagging the reader's 
mind, is whether one can actually perform our experi- 
ment. The answer is threefold. To begin with, nobody 
knows of the existence of a material exhibiting an SU (2)- 
order parameter structure. Fortunately, the existence 
of two spin superfluids is well-established: '^He-A and 
■^He-B. We will show that '^He-B has an order parame- 
ter structure similar to that of the pure spin superfluid. 
The effect of dipolar locking will destroy the spin vortex 
caused by the electric field, however, see Section (jXIBI) . 
Then we will show that "^He-A has, for subtle reasons, 
the wrong topology to perform our experiment. We will 
also demonstrate that the small spin-orbit coupling con- 
stant forces us to use an amount of ^He with which one 
can cover Alaska, turning our experiment into a joke. In 
the outlook of this work, we will discuss how the organic 
superconductors2^j2§ might meet the desired conditions. 
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Let us first consider the secrets of ^He more generally. 

XI. AND ORDER PARAMETER 

STRUCTURE 

As is well-known, '^He is a fermionic atom carrying spin 
i. In field theory, we describe it with an operator Cpa, 
where p is momentum and a is spin. In the normal phase, 
it is a Fermi liquid, but for low temperatures and/or high 
pressures, the He displays a BCS-like instability towards 
pairing. Indeed, the condensate wave function displays 
an order parameter which transforms under both spin 
and orbital angular momentum: 

(*| ^ pCp„C_p0 |*> = ^^.(iCT'^fT^)^^, (99) 

p 

so the order parameter describes a p-wave state. The 
A^i carry a spatial index i and an internal spin index ji. 
The numbers A^i transform as a vector under the spin 
rotation group S0{2>)^ acting on the index and the 
orbital rotation group iS'0(3)'^ acting on the index i. We 
can reconstruct the wave function \^ > from the A^j^i as 
follows. First we rewrite them as a vector decomposition 
with amplitudes a^i in the following way: 

Af,^ = Y,aki\i\[- (100) 

k,l 

The A'^'' are vectors. Then the wave function in momen- 
tum space 5'(p) =< pj^" > is the decomposition 

*(p) =^aHK=i,fc(p) Xs=M , (101) 

k,l 

where Yl^i k 

is a triplet spherical harmonic and Xs=i.z 
is a triplet spinor. This means that the order parameter 
has 3x3x2 real degrees of freedom. Indeed, follow- 
ing Volovik^° and Leggetlfii, there exist two mean-field 
states. 

The first one is an isotropic state with vanishing total 
angular momentum J = L -|- S" = 0. In order to have zero 
projection of the total spin nij = mi + rus = 0, we have 
for the coefficients in the decomposition (|100p 

a+- = a-+ = aoo = A^. (102) 

This state is called the B-phase of '^He, or the BW-state, 
after Balian and Werthamer— . This means that the or- 
der parameter looks like 

A^, = AgSar. (103) 

There is still a degeneracy, however. Indeed, both the 
spin and orbit index transform under SO(3), which leads 
to an order parameter manifold 

Rc.^^ Rf,R^aSa^, Or i? = i?^ (i?^ ) " ^ . (104) 



So the matrix R G SO{3) labels all degenerate vacua, and 
describes a relative rotation of spin and orbital degrees 
of freedom. Including also the U{1) phase of the matter 
field, the order parameter manifold of '^He-B is 

Gb - S0{3)rel X U{1) 

■matter • 

(105) 

This will be the starting point of our considerations for 
■^He-B, in which we will often drop the U{1) matter field. 

The second one is the A-phase, which has just one 
non- vanishing amplitude in (jlOOp . 

ao+ = V2Aa, (106) 

which corresponds to a state with rris = and mi — 1. 
The quantisation axes are chosen along the i-axis, but 
this is just arbitrary. This is known as the '^He-A phase, 
or the Anderson-Brinkman-Morel (ABM) state^l. The 
order parameter is 

Aai = AAZa{xi + iyi). (107) 

Rotations of the quantisation axis of ^He-A lead to the 
same vacuum, which tells us how to describe the degen- 
eracy manifold. The vector describing spin, called the 
d-vcctor in the literature''^-, can be any rotation of the 
z-axis: 

d„ = (108) 

Since only the direction counts in which the d-vector 
points, its order parameter manifold is the 2-sphere S^. 
The orbital part of the order parameter is called the I vec- 
tor, which is in the "gauge" Eq. ()107p simply z. Again, 
the orientation is arbitrary, so that any rotation R^ and 
gauge transformation e*"^ leads to a correct vacuum state, 

eW+*ef)=e'^4(£,+iy,), (109) 

where I = e^^-* x e^^-* is invariant under e**^. This 
phase communicates with the phase of the matter field, 
so that the order parameter has a relative U{l)rei — 
U{\-) matter -orbital- For the determination of the order 
parameter manifold for He-A, we need to observe that 
the order parameter does not change if we perform the 
combined transformation d —>■ —d and (e^'' -I- «e|^^) — *■ 

— (ep'' + iep'). This means that we have to divide out 
an extra Z2 degree of freedom. In summary, the order 
parameter manifold for He-A is 

Ga = [S^, X SOi3)i)/Z2, (110) 

where s refers to the spin and I to the orbit. The in- 
tricateness of the order parameter already indicates that 
there is a lot of room for various kinds of topological 
excitations and other interesting physics. For exten- 
sive discussions, we recommend the books of Grigory 
Voloviki2i2^. What counts for us, however, is how the 
topology is influenced by switching on fixed frame gauge 
fields. 
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^He-B 



B. Dipolar locking 



As discussed above, the order parameter of ^He is de- 
scribed by an 5*0(3) matrix R. The question is now if 
R admits spin vortex solutions. In principle, it does, 
because SU{2) rotations are like 5*0(3) rotations, since 
they are both representations of angular momentum, as 
we learned in freshman quantum mechanics courses. This 
means that, in principle, all considerations for the SU (2) 
case apply to ^He-B as well. In particular, the spin su- 
perfluid velocity Eq. (I75p has a similar expression, but 
now with g = R E 50(3). It reads 



-eapjRi3jd,R^j. 



(Ill) 



Inspired by the SU{2) case, which was effectively 
Abelianised, we try a vortex solution around the z-axis 
(assuming the electric field is radial) 



R = exp(i0J3) 




In the 1970s, Leggett described in his seminal article 
about '^He many important properties of this interest- 
ing systemii. One of them is how the spin part of the 
condensate wave function '^{x) interacts with its orbital 
motion by a 5 • £ interaction. According to Leggett, the 
contribution of the Cooper pairs to the dipolar energy is 

Ed^p = ~gd^pJ dx^ (|*(2^)|' " 3|f • *(f)|2) 

— 3|n • {Aaina)]"^ - constant, (115) 

remembering that the spin order parameters carry a spa- 
tial index, cf. Eq.'spI37)l. piI5)l . We used the notation 
n = On inserting the order parameters Ea.'s (|107p 
and (|103p . we obtain for both phases the dipolc locking 
Lagrangians 



(112) 



where J is the generator of total angular momentum, and 
e = arctan(l^). With the help of the 50(3) analogue 
of Eq. ([75|) . the superfluid velocities Ea. (|llip are readily 
calculated to be 



Ldip,B 
Ldip,A 



-gd^p{{^TRf + ^r{Rf), 
-gdip{l-df- 



(116) 
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— -{d2Rlk)R2k — — ~^^Ei 



-{d3R2k)Rik = 
-id3Rik)R3k = 0, 



(113) 



For the '^He-A part, we do not need to solve the equa- 
tions of motion in order to infer that the orbital and 
spin vector wish to be aligned. For the B-phase, we 
give a derivation of the Leggett angle. A general ma- 
trix R G 50(3) can be described by three Euler angles. 
For the trace, only one of them is important, let's say it 
is called 6. Then 

Ld^p,B = ~9d^p {(1 + 2cos6l)2 + 2(cos2 - sin^ 0)} , 



xj. Since the groups 50(3) and SU{2) which leads to the static equation of motion 

„ dLdip^B , n. , 

O = ^^=4cos0-1, 
with the Leggett angle as solution, 

0L = arccos(-i) ~ 104°. 
4 



where = : 

give the same equations of motion Eq. (|83p . we see that 
the Ansatz Eq. l|112p satisfies these as well, giving a spin- 
Hall current for the 0-polarised spin. In other words, 
in ^He-B is a possible candidate for our quantised spin 
vortex. 

This result can also be understood by topological 
means, in the following way. The equation of motion 
for the SU{2) case tells us, that the vacuum manifold 
for the spin becomes U{1) instead of 50(3) ~ SU{2). 
Only if we were allowed to change the orientation of the 
wire, described by a point on 5^, we would obtain the full 
50(3). This is the translation of the mathematical fact 
that 50(3)/5^ ~ U{1), merely saying that a rotation 
is fixed by an axis of rotation and the angle of rotation 
about that particular axis. The implication is that we 
need to calculate the fundamental group of Gb/S^ in- 
stead of Gb itself: 



(117) 



(118) 



(119) 



The Leggett angle tells us that one degree of freedom 
is removed from the order parameter of ^He-B so that 



50(3) 



rel 



G 



B,a 



(120) 



but 7ri(5^) — 0, as any closed path on the sphere can be 
continuously shrunk to a point. 

Now we can also understand that dipolar locking de- 
stroys vortices, even in a cylindrical set up, i.e. with a 
radial electric field, since 



7ri(50(3)/52) =7ri(C/(l)) = Z, 



(114) 



leading to the existence of vortices in a cylindrical set up, 
i.e., the inclusion of radial electric fields induces vortices. 

There is however one effect which destroys our spin 
vortex solution. This effect, known as dipolar locking, 
will be discussed in the next section. 



T^i{GB,dipl S'^) = 7ri(e) = 0. 



(121) 



The "division" by the manifold 5^ translates the fact 
that different vortices in the '^He-B manifold are only 
equivalent to each other up to different orientations of the 
cylindrical wire, being described by 5^. Another way to 
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FIG. 5: The destruction of the spin vortex by dipolar locking. 
The f/(l) degree of freedom is indicated by an arrow. In the 
center where the electric field is located, the angle follows a 
vortex configuration of unitwinding number, corresponding to 
one charge quantum. Since the electric field, decaying as i , is 
not able to compete with the dipolar locking at long distances, 
the f/(l) angle becomes fixed at the Leggett angle, indicated 
by a horizontal arrow. 



understand the destruction of vortices beyond the dipolar 
length, is that the U{1) vortex angle 9 is fixed to the 
Leggett angle, as depicted in figure O 

The fact that the vortices are destroyed, even though 
the spin-orbit coupling energy is higher than the dipolar 
locking energy^^, is due to the fact that small energy 
scales do play a role at large distances. This is similar 
to spontaneous symmetry breaking in, for example, an 
XY-antiferromagnet. A small external field is enough to 
stabilize domain walls at long wavelengths. 



C. ^He-A 

In the discussion of the pure spin superfluids and of 
•^He-B, we used the fact that the order parameter has 
a matrix structure, namely SU{2) and 5*0(3), respec- 
tively. For the SU (2) case we had to transform from the 
fundamental spinor representation to the adjoint matrix 
representation. Since both representations are SU{2), 
the physics did not change fundamentally. The resulting 
equations of motion were equations for group elements g, 
with the ramification that spin vortex states lower the en- 
ergy with respect to the trivial solution, cf. Eq. (PT|l . Asa 
result, the vacuum manifolds in both cases become U{1) 
instead of SU{2) (pure spin superfluid) or 5*0(3) (^He-B 
without dipolar locking). The topological protection of 
the spin vortex solution followed from the fact that C/(l) 
is characterised by the winding numbers, 7ri(C/(l)) = Z. 

For the case of '^He-A, matters are different, since the 
spin order parameter for ■^He-A is a vector in 5^ instead 
of a matrix in 50(3). Although 50(3) acts on 5^, these 
manifolds are not the same. What we will prove is that 
as a result, spin vortices do not lower the energy in the 
presence of an electric field, as opposed to the '^He-B and 
pure spin superfluids. The consequence is that the vac- 
uum manifold remains S'^, and since 7ri(5^) = 0, spin 
vortices are not protected. The presence of dipolar lock- 
ing will not change matters. 

Let us prove our assertions by deriving the equations of 



motion from the Lagrangian for "^He-A. The free energy 
functional^?- for ^He-A is quite analogous to that of a 
liquid crystal^, as the A phase is both a superfluid and a 
liquid crystal in some sense. Besides the bulk superfluid 
energy, there are also gradient energies present in the 
free energy, of which the admissible terms are dictated 
by symmetry: 

Fgrad = Jl{diAaj){diAaj)*+j2{diAat){djAaj)* 
+j3{d,Aa,){djA 

at) (122) 

= A^J„e'^'-(er^+«ef^). 
This then leads to 
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Fgrad °" — -^Kijmndiemdjen + Cij {Vs)iejkldkei 

+ ^Pi3{9ida){djda) + gdipidaEa)"^ . (123) 

The coefficients Kijmn and Cij are the liquid crystal like 
parameters-*^. The superfluid velocity Vs is the Abelian 
superfluid velocity coming from the relative U{1) phase. 

We are going to prove that "^He-A does not have topo- 
logically stable spin vortices, and that dipolar locking 
does not stabilize these. Generically, the spin stiffness 
tensor pij is given bj^ 



p 



(124) 



but it becomes fully diagonal when we neglect 
anisotropics in the spin wave velocities, i.e., p" = p-^ . 
We also assume that the Kij^mn and and Cij are fully 
diagonal, since this will not change the nature of the 
universal low energy physics. Including now spin-orbit 
coupling and kinetic terms the '^He-A Lagrangian is 



L^(i^„,,^,B) = - 



2mc^ 
2fimn 



[\doej\^ + {dod^f + 



2/iTOns 



^a/Bf^aikdpdidjEk 



l{E^-B^)-\g,.,{d.iy. 

(125) 

The strategy for solving the equations of motion is as 
follows: flrst we demonstrate that a spin vortex is pos- 
sible without dipolar locking, but that it does not gain 
energy with respect to the constant solution. Then we 
show that the spin vortex is not stabilized by switching 
on the dipolar locking. 

Without dipolar locking a spin-only action is ob- 
tained, leading to an equation of motion which resembles 
Eq.i 



ft 



did 



2fimns 



{ea(3j)di3Ek 



= 0. 



(126) 
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Let us choose a reference vector D,^, such that dj = 
Rj^D^. Again, R is an 5*0(3) matrix, describing the 
superfluid phase of the S'^ variable d. In this way, the 
equation of motion for the group element R reads 



diRjv 



2fimns 



ik{^a(3j)Rl3vEk 



0. (127) 



Using cylindrical coordinates, the demonstration that the 
spin vortex Ansatz for i? is a solution to this equation 
of motion is analogous to the proof that a spin vortex 
exists in ^He-B, cf. Ea. (|113p . On the other hand, this 
equation also admits a constant i?, i.e., Eq. p26p admits 
a constant Z)^ as well. Substituting both solutions back 
into the energy functional Eq. p25p , no energy differences 
between the spin vortex and the constant solution show 
up. In mathematical terms, the vacuum manifold in the 
presence of a cylindrical electric field remains S*^ . In plain 
physics language: the electric field does not prevent phase 
slips to occur. 

The presence of dipolar locking makes matters even 
worse, since the equations of motion become equations 
of motion for e and d involving dipolar locking, 

2/imns 



didj 



(128) 



It is clear that in general, a vortex configuration for d is 
not a solution, since the left hand side of the equation 
for d is annihilated, whereas the right hand side is not. 
Instead, the orbital and spin vectors will perform some 
complicated dance, set in motion by the electric field. 

The verdict: our charge trapping experiment will not 
work employing '^He-A. 



D. Baked Alaska 

In the search for an experimental realisation of the 
proposed charge trapping experiment, it turned out that 
■^He-B admits spin vortex solutions only at short wave- 
lengths. But if there were a way to circumvent dipolar 
locking in some ideal world, nothing would stop us from 
performing the actual experiment. 

Or... does it? It turns out that the numbers which 
Nature gave us, conspire to obstruct matters. It is really 
hidden in the fact that electric fields are so strong, and 
spin-orbit coupling so weak. Let us first confess that 
in the previous considerations, we did not regard a very 
important part of our charge trapping device, namely, 
the wire itself. The charge stored on it is hugely repelling 
indeed, giving rise to an enormous charging energy. 



First, we calculate the Coulomb energy stored in the 
wire. Let p{x) be the charge density distribution, which 
we approximate by a step function of the radius. Then, 



rrr 1 f Pi^) Pi^) , , I 

Wcouiomb = / \,, dxrfx' 



87reo J ||x — x'l 



1 Qlot 



STrep iTa?L 



(129) 



We integrated over the center-of-mass coordinate, and 
(with the definitions u = x — x' and r = L/a) we intro- 
duced 



duz I 2ndu±u±- 



1 



^2.l^-lL^^a£du.f7{^'^ 

= 27r|-iL2 + ^(^^yrT^ + ln(g+v/rT^))| 



27r— ln{2q) for L » a. 



(130) 



We used the standard integral / drVl + 
irVl + t2 + i ln(r + Vl + t^). Hence 
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W'coulomb = Z ^ -^1^ 

ovreo \ a 



(131) 



For the parameters under estimation, Wcouiomb/-^ — 
IJ/m, which is really enormous, since the coupling con- 
stant of electric fields is so huge. 

The question is now if the superfluid is strong enough 
to keep the charge trapped. Indeed, if it doesn't, the 
system can lower its energy by simply discharging the 
wire, causing a big spark, and destroying the superfluid. 
This is analogous to magnetic flux trapping in supercon- 
ducting rings with the Aharonov-Bohm effect'^^ . The flux 
trapped in a ring is a metastable state, but the supercon- 
ducting condensate is strong enough to keep it there. 

However, spin-orbit coupling is too weak to do so with 
our Aharonov-Casher analogue. In fact, the only thing 
the system needs to do, is to destroy the spin superfluid, 
not in the whole container, but just a small strip of the 
order of the coherence length ^, which is of the order of 
0.01/im'^^. We now need to estimate the energy density 
of the fluid. To do this, we perform Landau theory for 
the superfluid order parameter 



6F 



(132) 



This expression is zero when there is no superfluid. There 
is no kinetic term, since ip is parallel transported by the 
electric field: indeed, if it satisfies the equations of mo- 
tion, the kinetic term vanishes, cf. Eq. ([55)1 . Hence, 
we are only left with the potential energy terms. From 
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enough to cover Alaska, and much more than the total 
amount of He on Earth (180 liters). There might be 
enough He on the Moon, but still it is a "only in your 
wildest dreams" experiment. Is there no way out? In 
the concluding section, we give a direction which might 
provide some hope. 



FIG. 6: View from the top of our container. The container 
radius is R, and the wire has radius a. Now, the Coulomb 
energy of the wire has to make a tiny region of superfluid 
normal again, in order to make phase slips happen, removing 
the topological constraint. The region in which this should 
happen, needs to be of the width of the coherence length ^, 
but it has to extend over the whole radius of the container. 



Landau theory, we know the saddle point value for in 
terms of a = a(T — Tc) and 6, viz., 

m' = ^^6F = -V^{T-T,), (133) 

where V — ttR^L is the volume of the container. Note 
that R is the unknown variable in our problem. From 
Landau and Lifschitz we obtain the BCS-parameters 

ksTc 



67r2 kgT T\ 

a{T) = —— (fcsTJ 1-— , b^a 

7C(3) M \ TcJ p 

(134) 

where p is the superfluid density. For low temperatures 
T << Tc we have p ~ ep, s 



5F ~ 3.52{nkBTcV) 



(135) 



We use experimental values^ ep/kB = 0.312_ft' and Tc ~ 
3mK. From the Fermi gas relation p = p^/Sw^h^ we 
then obtain p « 15 mol/litcr. This leaves us with an 
estimate 



~V 



34 J/m^ 



The question we need to ask is: how big does the con- 
tainer radius i? need to be, in order to remain in the 
mctastable, charge trapped state? Per length unit L, the 
estimate is 



^Coulomb SF 



(136) 



Due to the enormously small ^ and the enormously big 
Wcouiomb, this leads to a truly disappointing radius of 



R ~ lOOO/cTO, 



(137) 



XII. OUTLOOK: ORGANIC 
SUPERCONDUCTORS 



In the previous section, we have seen that the small 
spin-orbit coupling energy and the big electric fields are 
disastrous. This is due to the fact that the coherence 
length ^ is small. In turn, the reason for that is that in 
Landau theory, ^ oc In other words, the heavier the 
constituent particles, the worse things get. So we need to 
look for lighter things. The first candidate would be elec- 
trons, since they are 5000 times lighter. However, as they 
are charged, charge effects highly overwhelm the whimpy 
spin-orbit coupling effects. So we need something made 
out of electrons, having however a huge gap for charge 
excitations: we need a spin superfluid made out of a Mott 
insulator. Does this exist? 

In recent years, there have been many advances in 
the research on highly frustrated systems on triangular 
lattices^, which are realised in organic compounds. In 
the last two years, Kanoda et al. have done specific heat 
measurements in the spin liquid phase of the organic su- 
perconductor K-(ET)2Cu2(CN)3, see Figure[71 Although 
the spin liquid state is known to be a featureless para- 
magnet, the specific heat shows a linear behaviour as a 
function of temperature^ii^^. 

The linear behaviour has led theorist P.A. Lee to the 
idea that this might be caused by fermionic spinons form- 
ing a Fermi surface'^''. It is plausible that at low energy 
scales, a BCS-like instability of the Fermi surface might 
give rise to an S* = 1 spinon condensate. This would 
then be the desired spin superfluid made out of a Mott 
insulator. The theoretical complication is that due to the 
SU{2) slave theories developed by Lee and Wen^^, there 
will be transversal gauge degrees of freedom, blocking 
the triplet channel, which should give rise to some scep- 
ticism about whether the organics are able to become a 
triplet superfluid. Whether or not this is the case, to our 
opinion, the idea of charge trapping provides a good mo- 
tivation to pursue the BCS-instability towards a triplet 
state of the spinon metal further. 
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FIG. 7: The phase diagram of the highly frustrated k- 
(ET)2Cu2(CN)3, as proposed by Kanoda^^ The spin liquid 
state shows linear specific heat, which might signal the pres- 
ence of a spinon Fermi surface. This would amount to making 
a spinon Fermi liquid out of an insulator. Then the interesting 
possibility is that this spinon metal might be unstable against 
an 5* = 1 spin superfluid. 



and substitute them in the Pauli equation Eq.([9]). In or- 
der to obtain Eq. ([25|) . we multiply the Pauli- like equation 
by r'^/2 and use the expressions 



APPENDIX A: USEFUL FORMULAS 



In order to obtain Eq's(l2T|) and ([26]) we calculate the 
expressions 
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